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Abstract 

An analysis of inelastic photoproduction of J/ip mesons is presented using data collected 
at the ep collider HERA corresponding to an integrated luminosity of above 80pb _1 . Dif- 
ferential and double differential cross sections are measured in a wide kinematic region: 
60 < W yp < 260 GeV, 1 < p\ ^ < 60 GeV 2 and 0.05 < z < 0.9, where z is the fraction 
of the energy of the exchanged photon transferred to the J / ip meson in the rest frame of the 
target proton. Cross sections at z $ 0.3 are presented for the first time. Theoretical calcu- 
lations within the Colour Singlet Model at NLO for direct photon processes are shown to 
give a good description of the data in the medium z region (0.3 < z < 0.9) up to the high- 
est pf^ values. A calculation using a k t factorisation approach in LO in the Colour Singlet 
Model is also able to describe these data. The data in the full z range are also compared 
to LO calculations within a non-relativistic QCD framework including colour octet and 
colour singlet contributions for direct and resolved photons. It seems possible to reconcile 
data and theory with modest contributions from colour octet processes. The polarisation of 
the J/ij) meson is measured as a function of z and pt^ and is reasonably described by the 
theoretical predictions. 
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1 Introduction 



The production of J ftp mesons, ep — > eJ/ipX, has been studied intensively at HERA [[T|-|4|]. 
The high available energy allows the contributing mechanisms to be studied in a wide kinematic 
range in both Q 2 and W 1P . Here Q 2 is the negative squared four-momentum of the exchanged 
photon and W 1P is the centre of mass energy of the photon proton system. In this paper we 
concentrate on the inelastic process where, in the proton rest system, the J/ip meson carries a 
fraction z < 0.9 of the photon energy in contrast to diffractive processes where this fraction is 
close to 1. The inelastic process is dominated by boson gluon fusion: a photon emitted from the 
incoming lepton interacts with a gluon from the proton to produce a cc pair which subsequently 
forms the J/ip meson. Here we will consider only interactions of quasi-real photons|], Q 2 ~ 0. 
In this case the photon can either couple to the c quark directly ("direct" processes, Fig. [TJa or 
b) or can interact via its hadronic component ("resolved" processes, Fig. [TJc). 



a) e e b) e e c) e 




Figure 1: Generic Feynman diagrams for inelastic J ftp production. a,b) direct photon pro- 
cesses; c) resolved photon process. In diagrams a) and c) the cc pair leading to the formation 
of the J ftp can be in a colour singlet or octet state while in b) it can only be in a colour octet 
state. Additional soft gluons emitted during the hadronisation process are not shown. 



Many models have been suggested to describe inelastic J ftp production in the framework of 
perturbative Quantum Chromodynamics (pQCD) such as the Colour Singlet Model (CSM) 0], 
the colour evaporation model [Q] and soft colour interactions [|7p. Recently, discussions have 
focussed on calculations based on a factorisation theorem in a non-relativistic QCD approach 
(NRQCD) [Bp. These theoretical descriptions differ in many details. One important difference 
concerns the states they allow for the intermediate cc pair produced in the hard boson gluon 
interaction in terms of angular momentum and colour. Furthermore they differ in the non- 
perturbative description of the transition from the intermediate cc pair to the J/ip meson. 

In the CSM it is assumed that the intermediate cc pair is only produced in the quantum state 
of the J ftp meson, i.e. in a colour singlet state with spin 1 and no orbital angular momentum. 
However, the CSM predicts prompt production rates of J/ip and ip' mesons in pp interactions, 
which are lower by more than an order of magnitude than the rates seen by the CDF collabora- 
tion [0] . The NRQCD approach, in which contributions from cc pairs in colour octet states are 
predicted in addition to colour singlet contributions, can accommodate the measured pp cross 
sections. Long distance matrix elements (LDMEs) are used to describe the transition of the 
various intermediate cc states to the J/ip meson. In previous publications []2|,|3|] it was shown 

'The regime of virtual photons is studied in [ID]. 
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that in photoproduction at HERA the colour octet contributions are not observed with the mag- 
nitude predicted by leading order NRQCD calculations, in which LDMEs are used which are 
extracted, also in leading order (LO), from pp data. However the data were well described by 
the next-to-leading order (NLO) calculations in the CSM. Subsequently theoretical efforts were 
made to improve the simultaneous description of the photoproduction and pp processes within 
NRQCD, e.g. the effects of higher order corrections have been estimated in the extraction of 
the LDMEs as well as in the photoproduction amplitudes jrDj— |T~3|] . 

Since our first publication [|2|] the amount of data taken by the HI experiment has increased 
and allows more detailed investigations in a larger region of phase space. Here we present 
an analysis in an extended z region, 0.05 < z < 0.9 covering 60 < W 1P < 260 GeV and 
transverse momenta of the J/ip, 1 < p\^ < 60 GeV 2 . In the medium z region, which has 
been analysed previously, mainly direct photon processes are expected. At the lowest z values 
resolved processes are expected to contribute to the cross section in addition. 

The data in the medium z range are compared with two sets of calculations within the CSM, 
namely a NLO calculation [ |T?| ] and a LO Monte Carlo computation [ |T51 , |T6T ] using a "k t factori- 



sation" approach^ [|T7[|r8|] . A comparison of the data with LO calculations in the NRQCD frame- 
work including colour octet contributions is carried out over the full range 0.05 < z < 0.9. The 
decay angular distribution of the J / ip meson is analysed yielding information on the polarisa- 
tion of the J/ip meson and the results are compared with predictions in the NRQCD framework 
both with and without colour octet contributions [ |19| ] and in a k t factorisation approach 



2 Theoretical Models for J /ifr Production 

The salient features of the relevant theoretical calculations are briefly summarised here. A 
comprehensive overview can be found in [[H]] . 

In the Colour Singlet Model [|]] for direct photon processes the J/ip meson is produced via 
75> — > cc [l, 3 5i] + g, where the quantum state of the cc pair is described in spectroscopic nota- 
tion!]. The transition from the cc pair to the J / ip meson is calculated in potential models and can 
be related to the measured leptonic decay width of the J / ip meson. For photoproduction, full 
NLO calculations, i.e. up to order a^, where a s denotes the strong coupling, are available for the 
direct process [|14|]. The NLO contributions lead to an increase of the cross sections by roughly 
a factor 2, dependent on the transverse momentum of the J/ip meson. The resolved photon 
process is expected to be dominated by gg — > cc [l, 3 <Si] + g, for which only LO calculations 
are available. 

These CSM calculations are performed using standard collinear gluon density distributions. 
An alternative k t factorisation approach [JTV|, P~8l] has recently been successfully applied to the 
description of a variety of processes [fT5|]. In this approach the J/ip production process is fac- 
torised into a k t dependent gluon density and a matrix element for off-shell partons. A LO 



calculation within this approach is implemented in the Monte Carlo generator CASCADE [|T3|]. 



2 Here "fc t " refers to a transverse momentum component of the gluon entering the hard interaction. 

3 Spectroscopic notation: 2S+1 Lj where S, L and J denote the spin, orbital and total angular momenta of the 
cc system that is produced in the hard process. The first number in the square bracket indicates the colour state of 
the cc pair. 
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In the NRQCD approach the J /if) production amplitude factorises into short distance terms for 
the creation of colour octet or singlet cc states with definite angular momenta and long distance 
terms describing the transition of these cc states to the J /if) meson. A double expansion in the 
strong coupling parameter a s and v, the relative velocity of quark and antiquark, is obtained. 
In specific calculations [|TD|, [TT], ^T]] only the most important contributions are kept. In this 
expansion the lowest order term in v is the colour singlet term. Assuming that all other terms 
do not contribute, the CSM is recovered. 

Whereas the short distance amplitudes are calculable in pQCD the LDMEs must be determined 
from experimental data at present. They are assumed to be process independent, i.e. they can 
be determined for example in pp collisions or from B decays and can then be used in predic- 
tions for photoproduction. Calculations for photoproduction based on the NRQCD approach are 
available in LO taking into account the contributions [l, 3 Si], [8, 3 Si], [S^So] and [8, 3 Pj =0 ,i,2]- 
The relative strength of the colour octet contributions depends crucially on the size of the corre- 
sponding LDMEs. Unfortunately the values for the LDMEs important at HERA still show large 
uncertainties (a summary can be found in [|2~T]]). Early extractions of the LDMEs from pp data 
were performed in LO. Subsequent estimates of higher order effects led to considerably smaller 
values [0-jnp. 

The first comparisons of NRQCD calculations with photoproduction data [0,|3|] revealed large 
discrepancies between data and theory at high values of z, the theory being a factor 3-5 higher 
than the data. In later calculations resumming of soft gluon emissions within the NRQCD 
expansion [ [TT| ] was shown to damp the colour octet contributions at high z. The region of 
validity of these calculations is restricted to J /if) production with high transverse momentum, a 
region where experimental data suffer from limited statistics. 

All theoretical calculations contain a number of free parameters, e.g. the parton density distribu- 
tions, the values of a s and the charm quark mass m c as well as the choice of the renormalisation 
and factorisation scales. In the NRQCD approach the values of the octet LDMEs are additional 
parameters. The comparison with the data in the NRQCD approach also suffers from the un- 
certainties associated with LO calculations. Although the NLO terms have not been calculated 
in the NRQCD approach similar effects as in the CSM may be expected, where the NLO terms 
lead to an increase of the cross section of typically a factor 2 with a strong p t ,^ dependence. 
Measurements of the polarisation of the J/ip meson should help to discriminate between differ- 
ent theoretical descriptions, independently of normalisation uncertainties. 

Table [I] gives an overview of the calculations available in the low and medium z regions. 



z region 


Contributing Process 


Available theoretical calculations 


medium z 
(0.3 <z< 0.9) 
low z 

(0.05 < z < 0.45) 


direct photon process 

direct + resolved photon processes 


CSM NLO and LO, NRQCD LO 

k t factorisation 

CSM LO, NRQCD LO 



Table 1 : Overview of theoretical calculations compared with the present data in the medium 
and low z regions. 
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3 Experimental Conditions 



The data presented here were collected in the years 1996-2000 and correspond to a total in- 
tegrated luminosity of (87.5 ± 1.3) pb _1 . HERA was operated with electrons or positrons of 
27.5 GeV. The proton beam energy was 820 GeV before 1998 and 920 GeV since. 



3.1 Detector and Event Selection 

The experimental methods are similar to those described in [Jl|, H»^2[]- Details of the detector 
and the analysis can be found in [ |2~3"| ] and [E3J], respectively. J/tp mesons are detected via the 
decays J/tp — > (branching fraction of (5.88 ± 0.10)% [|26*|]). Reconstructed tracks in the 

central tracking detector (CTD) are identified as muons either because they produce minimal 
ionisation in the main (liquid argon, LAr) calorimeter or because they are linked to a track in 
the instrumented iron return yoke of the magnet (muon detector). The trigger requirements rely 
on signals from track chambers and the muon detector, so for consistency at least one muon has 
to be reconstructed in the muon detector. 

The event sample is subdivided into three (partially overlapping) datasets to facilitate specific 
analyses. The datasets differ in the kinematic requirements for z, W rp and the minimum muon 
momentum p™ m imposed in the selection (see table 0). Dataset I is a medium z selection, 
0.3 < z < 0.9. Dataset II emphasizes the low z region and the selection imposes a restriction in 
the polar angle of the muons to suppress background at z £ 0.2. The z range, 0.05 < z < 0.45, 
is chosen to overlap with dataset I. The W 1P ranges for the datasets I and II are adjusted to 
give high acceptance in the corresponding z range. This results in a reduced W 1V range for 
dataset II. In order to determine the differential cross section da/dz in the full z range, a third 
sample, dataset III, is defined, which covers the same W 1P region as dataset II but covers a 
complementary z range identical to that of dataset I. 

In selecting J/tp events two identified muons with momenta^ > 1.1 GeV are required (datasets 
I and III) and p u > 0.8 GeV for dataset II. One of the muons must be reconstructed in the muon 
detector and has to have a momentum p^ > 1.8 GeV. In addition to the decay muons at least 
three additional tracks are required in order to suppress diffractive J/tp contributions. Fur- 
thermore a cut on the transverse momentum of the muon pair, pt^> 1 GeV, is applied. The 
restriction to photoproduction is made by removing events with a cluster of energy above 8 GeV 
in the electromagnetic part of the calorimeter. The mean Q 2 is then 0.05 GeV 2 . 

The distribution of the invariant mass of the two muon candidates after these cuts is shown for 
datasets I and II in Fig. 0. There is a considerable non-resonant background in regions of low 



Dataset 


z 


W JP [ GeV] 




p min [ Ge yj 


I 


0.3 - 0.9 


60 - 240 


20-160 


1.1 


II 


0.05 - 0.45 


120-260 


20 - 140 


0.8 


III 


0.3 - 0.9 


120-260 


20-160 


1.1 



Table 2: Overview of kinematic ranges and selection criteria for the three datasets. All datasets 
are restricted to Q 2 < lGeV 2 andp t ^ > lGeV. 
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2.5 3 3.5 4 2.5 3 3.5 4 

[GeV] [GeV] 

Figure 2: Distribution of the invariant mass of the muon pairs after all selection cuts a) in 
the range 0.3 < z < 0.9, 60 < W 7P < 240 GeV and pt^ > 1 GeV (dataset I) and b) in 
the range 0.05 < z < 0.45, 120 < W 1P < 260 GeV and p t ^ > 1 GeV (dataset II). The data 
points with statistical error bars show the unlike sign combinations, the histograms the like sign 
combinations. The lines are fits as described in the text. 

z due mainly to misidentified leptons as can be seen from the distribution of like sign pairs in 
Fig. 0b. The following procedure is adopted to extract the number of signal events in a given 
analysis bin: the mass distribution is fitted in the interval 2.45 < < 3.45 GeV with a 

superposition of a Gaussian to describe the signal and a function linear in mass to describe the 
background. The width and position of this Gaussian are determined from an overall fit to the 
data and are then fixed. The number of signal events is obtained by counting the number of 
muon pairs in the interval 2.9 < M w < 3.3 GeV and subtracting the fitted background in this 
interval. By varying the shape assumed for the background the systematic error in the number 
of signal events is estimated to be 5-15% at low z values and less than 1% at higher z values. 

3.2 Kinematics 

The following kinematic variables are used to describe J/ip production: the square of the ep 
centre of mass energy, s = (k + p) 2 , the negative squared four-momentum transfer Q 2 = 
—q 2 and the photon proton center of mass energy W 1P = yjp + q) 2 - Here k, p and q are 
the four-momenta of the incident lepton, proton and the exchanged photon, respectively. The 
scaled energy transfer y = q ■ p/k ■ p is the energy fraction transferred from the lepton to the 
hadronic final state in the proton rest frame. Neglecting the lepton and proton masses, a good 
approximation is W 2 ~ ys in the limit of Q 2 ~ considered here. The fraction of the photon 
energy transferred to J ftp meson in the proton rest frame is given by 

Pi> -P 
z = — — , 
q-p 

where p^ denotes the J/ip four-momentum. 



8 



The kinematic variables are reconstructed using the Jacquet-Blondel method p5|]: 



y = wX {E ' Pzh 

had 

(E-p z )j/Tp 
2yE e ' 

where E e denotes the energy of the incoming lepton. (E — p z )j/^ is calculated from the J/ip 
decay muons and had runs over all final state hadrons including the J/ij) meson. In the 
calculation of the final state hadronic energy a combination of tracks reconstructed in the CTD 
and energy depositions in the LAr and SpaCal calorimeters is used. 

The resolution for W 1V and p\^ is much smaller than the width of the chosen analysis intervals. 
The absolute resolution for z is typically 0.08 for z > 0.4 and 0.02 at z ~ 0.1. The size of the z 
bins is chosen so that in the Monte Carlo simulation the fraction of events reconstructed in the 
bin in which they were generated is above 60%. 



3.3 Comparison of Data and Simulation 

The Monte Carlo program EPJPSI [ |27| ] is used to generate two types of simulated datasets. The 
first contains J/ij) events from direct photon gluon fusion, the second from gluon gluon fusion 
corresponding to the main expected contribution to the resolved photon component. Both pro- 
cesses are calculated in the Colour Singlet Model in LO. Relativistic corrections are included, 
which enhance the cross section at high z values [E7p. The proton and photon gluon densities in 



the simulation are taken from MRS(A') [g8|] and GRV-LO [gg], respectively. The QCD renor- 
malisation and factorisation scales are chosen as Mj/^. Parton showers are simulated but their 
effect is small. 

The procedure to extract cross sections relies on the detailed simulation of the HI detector re- 
sponse. Particular care is taken to adjust the efficiency functions for the muon identification and 
trigger algorithms in order to achieve a reliable description as a function of the muon momenta 
and angles. This adjustment is performed using high statistics datasets from diffractive J/ij) 



production [24]. Remaining differences between data and simulation are taken into account in 



the systematic error. 

The two simulated datasets for direct and resolved photon processes are merged and normalised 
to the data, as follows (see also Fig. |]a-c and d-f for the medium and low z datasets I and 
II, respectively). First, the generated events of the direct photon simulation are reweighted in 
Pt^p according to a parametrisation of the measured distribution in the medium z range. The 
simulation is then normalised to the data in the range 0.45 < z < 0.9, where backgrounds 
from diffractive and resolved photon processes are expected to be negligible. Extrapolation of 
the simulation into the low z region shows that 35 J ftp mesons are expected in the interval 
0.05 < z < 0.15 whereas 80 J/ip meson candidates are observed. This excess of observed 
events over the direct photon simulation is attributed to the resolved photon component, which 
is normalised accordingly. In principle a contribution from bb pairs decaying to J / ip mesons 
should also be taken into account. However, the theoretical and experimental uncertainties are 
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still large for bb production, and the process is thus neglected here (see also discussion in the 
next section). 

A contribution of resolved photons of about 50% at 0.05 < z < 0.15 is obtained falling to 10% 
for 0.15 < z < 0.3 and being negligible for z > 0.3 (Fig. |3|e). Note that the resolved component 
is not reweighted in p\^ since good agreement of the superposition of the two simulations 
with the data is found (Fig. ^f). This may, however, be fortuitous because of the neglected 
b contribution, which is expected to be important particularly at high values of p\^. Current 
simulations (AROMA [^Up and EPJPSI [|Z7|]) normalised to the measured value for ep — > ebbX 
with Q 2 < 1 GeV 2 of (14.8 ± 1.3j^;|)nb 0, predict b contributions of the order of 50% or 
more in the highest analysis interval 20 <p 2 ^< 40 GeV 2 . 



3.4 Other Contributions to J/ ij) Production 



After the selection cuts described in section ^4] the J ftp sample is dominated by inelastic J/ip 
production, in which the J ftp is directly produced from the cc pair in the boson gluon fusion 
process. However, contributions remain from both the diffractive and inelastic production of 
ip(2S) mesons, the production of b flavoured hadrons and possibly also from \c states. The 
estimated amounts of these feeddown processes are summarised in Table |3|. Note that these 
contributions are not subtracted, either because they are negligible, such as diffractive back- 
grounds, or, because the dependence on the kinematic variables has not been measured and is 
only poorly known from theory, as in production of bb pairs with subsequent decays to J/ip 
mesons. 

Contributions from diffractively produced J ftp and ip(2S) mesons are effectively suppressed 
by the requirements of at least three tracks in addition to the decay muons and p 2 ^ > 1 GeV 2 . 
The remaining contribution from the cascade decay of ip(2S) is estimated with a simulation 



of diffractive ip{2S) production using the DIFFVM generator program Q32[]. DIFFVM is nor- 
malised to the observed signal of ip(2S) — > J/ipn + n~ in a sample of events with exactly four 
tracks. A residual contribution from diffractive 4>(2S) production of about 3% in the highest z 
bin is inferred. Averaged over the total z range this background is less than 1% in the medium 
z region and negligible at low z. 



Source 


Contribution (%) 




Dataset I, III 


Dataset II 


Diffractive ip(2S) -»• J/ip + X (0.75 < z < 0.9) 


1 (3) 




Inelastic ^(25) -> J/if> + X 


15 


15 


b ->• J/tj) + X (0.3 < z < 0.45) 


5 (19) 


25 


Xc -> J/ip + 1 


1 


7 



Table 3: Estimated contributions from other processes to the measured J/ip+X data sample 
in addition to inelastic boson gluon fusion. The amounts in parentheses refer to the z bins in 
parentheses. 
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The inelastic photoproduction of ip(2S) mesons with subsequent decays to J /if) mesons is ex- 
pected to contribute 15% to J/ip production [ |14| ] and to show similar dependences on the kine- 
matic variables because the production processes are the same. 

The contribution from the decay of b flavoured hadrons can only be estimated with large uncer- 



tainties, since few measurements exist [ plQ and these considerably exceed current theoretical 
estimates. Using the EPJPSI generator program and normalizing the production cross section 
to the measured value for ep — > ebbX with Q 2 < 1 GeV 2 [[3lj], a contribution of (25 ±6)% due 



to b — ► J/ip + X is estimated in the kinematic region of dataset II (see also Fig. |3|). At higher z 
and lower W 1V the background fraction is smaller. For example in the lowest z bin of dataset I 
a contribution of (19 ± 5)% is found which would correspond to 5 ± 1% in the total z region of 
dataset I. 

A further potential contribution comes from Xc decays, Xc j + J/ip. The Colour Singlet 
Model predicts that Xc mesons can only be produced in resolved photon processes due to quan- 
tum number conservation. From EPJPSI 7% of the total resolved contribution is estimated to 
arise from this source. The production of Xc mesons is possible also in direct photon gluon 
fusion via colour octet 3 Si states, which would result in a contribution of less than 1% in the 



medium z region [21]. 



3.5 Systematic Uncertainties 

The estimated systematic errors on the cross sections are listed in Table |]. The uncertainty from 
the track reconstruction is based on an error of 2% for the efficiency per track and takes into 
account the observed track multiplicity distribution. The systematic error of the total hadronic 
energy includes the uncertainty of the calorimeter energy scales (4% for the LAr, and 10% 
for the SpaCal) and leads to cross section error estimates of 4% (6%) in the medium (low) z 
regions. The systematic errors in the muon identification and trigger efficiencies are estimated 
from the residual differences between data and simulation. The error in the determination of the 
event numbers in the analysis bins is estimated by changing the description of the non-resonant 
background from a linear dependence on the invariant mass to a power law. The error on the 
acceptance due to restricting the polar angle of the decay muons is obtained by using different 
proton and photon parton densities which influence the W rp distribution. In general the changes 
in the acceptance are small. A significant uncertainty is found only at large W 1P in dataset I and 
at small W 1P in dataset II. The uncertainty in the acceptance at low z due to the modelling of 
the z dependence is found to be 4%. The low z dataset II contains a sizeable contribution from 
b — ► J/ip + X. The effect on the acceptance is estimated using the EPJPSI simulation of bb 
production with a relative contributions as in table The decay angular distribution of the 
J/ip meson has an influence on the acceptance mainly in the medium z region. Varying the 
decay angular distribution of the J / ip meson within the measured uncertainties^ changes the 
acceptance by ±8% nearly independently of other variables. 

4 The parameter A in equation ([!]) in section is varied from (used in the simulation) to +1 or —0.5. 



11 



Source 


Uncertainty [%] 




Dataset I, III 


Dataset II 


Track reconstruction efficiency 


6 


6 


Muon identification efficiency 


3 


7 


Calorimeter energy scales 


4 


6 


Signal events 


< 1 


5-15 


Trigger efficiency 


6 


10 


Acceptance: PDFs 


4-8 


5-10 


z dependence 




4 


J/ip decay ang. distr. 


8 


7 


b -> J/ip + X 


< 1 


3-12 


Integrated luminosity 


1.5 


Branching ratio 


1.7 


Total 


13.3-15.5 


18.9 - 24.3 



Table 4: Systematic uncertainties for the production cross sections in the medium z and low z 
regions. The error due to the J/ ip decay angular distribution does not apply to the polarisation 
analysis. 



4 Results 

Cross sections for ep — > eJ/ipX are converted to ^p cross sections using the equivalent photon 



approximation [ p3|p4| ]. The photon proton cross section is defined by: 

da ep = a w fy/ e (y)dy, 



where f y / e denotes the photon fluxP] integrated over Q 2 from the kinematic limit of Q 2 nin 
to the upper limit of the measurement, Q 2 max = 1 GeV 2 . The measured cross sections are 
given at the weighted mean for each bin assuming specific functional forms to fit the data, i.e. 

cr 7P (W 7P ) oc WjL and da/ dpi ^(x. (p 2 ^ + M|)" n as described below. 

The measured cross sections are displayed in Figs. |]-^]and the numerical values are collected 
in Tables |^-|TT|. In the next sections the measured cross sections are compared with theoretical 
calculations (Table [I]). 



4.1 Comparison of Cross sections with the Colour Singlet Model 

The photon proton cross sections a 7P , da/dz and do /dpi ^ in the medium z range are given in 
Figs. |] and |] and Table |[ The data are first compared with the NLO CSM calculations [|2l]]. 
Although the p 2 ^ range of the measurement has been considerably enlarged compared with 

5 ff/e(y) = £ (2m*y - + log , where Q 2 min = m\£- y , taken from [0. 
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earlier analyses the agreement between data and theory remains good (Fig. |]c). The band on 
the theoretical calculations in Figs. |] and |5] represents the uncertainties due to the variations of 
the mass of the charm quark, m c = 1.4 ± 0.1 GeV and of a s {M z ) = 0.1200 ± 0.0025 (see 
Table [7] for the other parameter values). In the W JP and z distributions of Fig. |]a and b these 
uncertainties affect mainly the normalisation. The shapes and normalisations of the calculated 
cross sections as functions of W 1P and z are in approximate agreement with the data for the 
two parameter sets, (ra c , a s ) = (1.3 GeV, 0.1175) and (1.4 GeV, 0.1225). The combinations 
m c = 1.3 GeV and a s (M z ) = 0.1225 (upper limit of the band) and m c = 1.5 GeV and 
a s (M z ) = 0.1175 (lower limit) can be excluded here. 

The NLO CSM calculation describes the p 2 ^ distribution rather well (Fig. |]c). This is not the 
case for the LO CSM calculation, which lies above the data at low p 2 ^ and falls too steeply 
towards higher values of p 2 ^. The p 2 ^ dependence is sensitive to the choice of m c and a s . The 
calculations with (m c , a s ) = (1.3 GeV, 0.1175) and (1.4 GeV, 0.1225), with which the W 1P 
and z distributions are well described, also give a reasonable description of the p 2 ^ dependence 
apart from a tendency to undershoot the data at high p 2 t ^ . The p 2 t ^ distribution in bins of z (Fig. [5| 
and Table ^ is also described well by the NLO CSM calculation. This is interesting because 
the size of the colour octet contributions in the NRQCD approach is predicted to depend on z. 

4.2 Parametrisations of Transverse Momentum Distributions 

The differential cross sections do /dp\^dz (Fig. |5] and Table |^) in the different z regions are 
found to be well described by a functional form (p 2 ^ + M|)~ n . Fits to the data have been 
performed and the results for the exponents n are summarised in Table |5|. There is a tendency 
for n to decrease towards lower z values. A similar analysis has been carried out for inelastic 
J /if) production in the range 2 < Q 2 < 100 GeV 2 (0.3 < z < 0.9, 50 < W 1P < 225 GeV) 
There the differential cross sections as a function of the transverse momentum in the photon 
proton center of mass system were fitted to the same functional form yielding a value n = 
4.15 ± 0.50, which agrees very well with the present result. 



W 1P [ GeV] 


z 


n 


60 - 240 


0.3-0.9 


4.49 ±0.15 


60 - 240 


0.75 - 0.90 
0.6-0.75 
0.3-0.6 


4.8 ±0.2 
4.6 ±0.2 
4.4 ±0.2 


120 - 260 


0.05 - 0.45 


4.1 ±0.2 



Table 5: Results of fits to differential cross sections da/dp 2 ^ of the form (p 2 ^ + M|) _n in dif- 
ferent kinematic regions. Total experimental errors (statistical and systematic added in quadra- 
ture) have been used in the fits. 

4.3 Comparison of Cross sections with NRQCD Calculations 

The comparison of the measured cross sections da/dp 2 ^ and da/dz with the calculations within 
the NRQCD approach are shown in Figs. ||-f| For this comparison data in the low z region are 
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included (Tables [R)|and |TTJ). There are considerable uncertainties in such LO NRQCD calcula- 
tions. The band in Figs. ^|and [7jgives an estimate from [ |2T| ] of the main uncertainty, which arises 
from the LDMEs obtained from J /if} production in pp interactions. For the upper limit of the 
band the LDMEs extracted in LO are used, while for the lower limit values including estimates 
of higher orders in the extraction are used [[T2]]. The distributions of p\^ are shown in Fig. 
and b for the low and medium z ranges, respectively. In both z regions the NRQCD calcula- 
tion pT] ] and data are compatible within experimental and theoretical uncertainties. However 
there seems to be a difference in shape: the data have a slightly harder spectrum than predicted. 
For the low z range one should keep in mind that the calculation is for charm only, while the 
data contain feeddown from b quarks, which is expected to contribute a harder transverse mo- 
mentum distribution. In both z regions the colour singlet contribution alone falls significantly 
faster than the data. 

The differential cross section da/dz extending over the full z range, 0.05 < z < 0.9, is com- 
pared with the same LO NRQCD calculation pT| ] in Fig. [7j for a cut p\^> 1 GeV 2 . Resolved 
contributions are found to dominate in the calculation below z £ 0.15 (0.3) depending on the 
choice of LDMEs. In the comparison of data and theory one should note that the contribution 
from b decays in the data is sizeable at low z (of the order of 25%). Nevertheless one may infer 
from Fig. [7J that the LO NRQCD calculation, including CS and CO contributions, is able to give 
a fair description of the data both in shape and in normalisation, if the LDMEs are chosen to 
be close to the lowest available estimates (lower edge of the shaded band). Using larger LDME 
values in this calculation leads to a strong increase of the theoretical cross sections at high 
and low z values which is inconsistent with the measurements. The colour singlet contribution 
alone, which is also shown separately in Fig. [7J is roughly 30% below the data for z £ 0.5 
although the shape is adequately described in this region. At lower z values the colour singlet 
contribution falls below the data by up to a factor 3. 

Although the uncertainties in the calculations are substantial, the measured cross section in 
the lowest z bin, 0.05 < z < 0.15, suggests that a resolved photon contribution is present: 
Correcting for b and inelastic ip(2S) feeddown using the estimates in Table ^| (section |3~7i] ) and 



neglecting their uncertainties the measured cross section is found to be 2 standard deviations 
above the direct photon contributions alone (including CS and CO contributions with a range of 
LDME values as in Table Q). 

A different calculation in the NRQCD framework has been carried out by Kniehl et al. [[H]] and 
is shown in Fig. [7J as 'HO improved'. In this calculation higher order effects were taken into 
account approximately using NLO parton density distributions for the photon and proton and 
using the LDMEs corrected for estimated higher orders effects (the parameters of the calculation 
are given in Table |7|). This calculation gives a good description of the shape of the data but 
a normalisation factor of 3, as suggested by the authors [JTO], is necessary to reconcile the 
predicted cross section with the data. 

The tendency for the colour octet contributions to rise towards high z values (upper end of 
shaded band in Fig. |7j) is at variance with the data as has also been noted previously This 
discrepancy may be due to phase space limitations at the upper limit of z where the emission 
of soft gluons in the conversion of the cc pairs to J /if} mesons is suppressed. This is not taken 
into account in [|TT)|,[2"T1]. In [ pT| ] a resummation of the non-relativistic expansion was carried 
out, leading to a decrease of the predicted cross section at high z. In Fig. |] (Table ||) the 
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measured cross sections da/dz for p t ^ > 2 GeV and for p t ^ > 3 GeV are compared with 
the results of these resummed calculations in the kinematic region of dataset I. The calculated 
curves have been roughly normalised to the data points at low z. The resummed calculation for 
A = 500 MeV gives an acceptable description of the data[] at pt^> 3 GeV while the agreement 
between data and calculation is still poor for p t ,^> 2 GeV or for lower A values. 



4.4 Comparison of Cross sections with k t Factorisation Calculations 

An entirely different approach to inelastic J ftp production within the CSM [fTo]] is implemented 



in the Monte Carlo program CASCADE [|T5]]. Here direct photon gluon fusion processes are 
computed in the k t factorisation approach [|T7|,|r8|] using an unintegrated (k t dependent) gluon 
density in the proton and the gluons may thus be off-shell. This gluon density was obtained from 
a fit [JT3J] to the HERA structure function data (using the CCFM parton evolution equations). 
Fig. [5] shows a comparison of the data with the results from CASCADE. In Fig. |9|a da/dz is 
shown as a function of z. Good agreement is observed between data and predictions for z % 0.8. 
At high z values the CASCADE calculation underestimates the cross section. This may be due 
to missing relativistic corrections which are not available for the off-shell matrix element. It 
could also indicate a possible missing CO contribution. The W rp distribution is shown for the 
restricted z range 0.3 < z < 0.8 (Fig. |9|b) and the CASCADE model is found to be in reasonable 
agreement with the data. An important improvement compared to collinear LO calculations is 
visible in the p\^ dependence (for 0.3 < z < 0.9, Fig. |^c), where the CASCADE predictions 
fit the data quite well. This is due to the transverse momentum of the gluons from the proton 
which contributes to the transverse momentum of the J/tp meson. 



4.5 Polarisation Measurement 

The polarisation of the J/ip meson is expected to differ in the various theoretical approaches 
discussed here and could in principle be used to distinguish between them, independently of 
normalisation uncertainties. The polarisation analysis is performed in the "target frame", the 
rest system of the J/ip meson, using the direction opposite to that of the incoming proton as 
reference axis z' . Two angles are defined: 9* is the polar angle of the positive decay muon with 
respect to the z' axis and $* is the angle between the plane of the decay muons and the plane 
defined by the photon and the z'-axis. 

The corrected data in the medium z range (dataset I) are normalised to the integrated cross 
sections and fitted to the forms flBfl: 



1 da 

a d cos 6* 



oc l + Acos 2 ^; (1) 



1# oc 1 + ^ + %os2$*. (2) 
3 3 V 

6 The parameter A describes the energy lost by the cc system, in its rest system, in the conversion to the J/tp 

meson. 
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z interval 


A 


V 


0.3 - 0.6 


1.1 ± 0.4 


-0.2 ± 0.5 


0.6 - 0.75 


0.6 ± 0.4 


-0.4 ± 0.5 


0.75 - 0.9 


0.1 ±0.4 


-0.6 ±0.4 


p t ^ interval [ GeV] 


A 


V 


1 - 2 


1.3 ±0.5 


-0.3 ±0.5 


2-3 


0.6 ±0.5 


-0.3 ±0.5 


> 3 


0.1 ±0.4 


-0.3 ±0.4 



Table 6: Fit results for the polarisation parameters A and u as functions of z and pt^. The 
errors are due to the total experimental uncertainties. 



The parameters A and v can be related to the polarisation of the J /if) meson. The cases A = 1 
and —1 correspond to transverse and longitudinal polarisation of the J /if) meson, respectively. 
Fits of the data to equations ([]]) and (Q) are performed in 3 bins of z or in 3 bins of p t ^. The 
polar angular distribution in z bins is shown as an example in Fig. [I7J. The results for A are 



listed in table g and plotted in Fig. [TTja and c. Positive values are preferred although a decrease 
is observed with increasing z and increasing p t ^. The fit of the $* distribution is also performed 
in z and p t ^ bins using the fitted values for A. The resulting v values are slightly negative (see 
Table |6| and Fig. |TT| b and d) and within errors independent of z and p t ^. 



In Fig. [IT] the results for three LO calculations are shown together with the data: the NRQCD 
prediction including colour octet and colour singlet contributions [JT9p, the colour singlet con- 
tribution alone and a calculation using a k t factorisation approach and off-shell gluons [|20|]. 
The errors in the data preclude any firm conclusions. None of the three calculations predicts a 
decrease of A with increasing z, while a decrease with increasing p t ^ as observed in the data 
is expected for the k t factorisation ansatz. The full NRQCD prediction is also compatible with 
this. In order to distinguish between full NRQCD and the colour singlet contribution alone, 
measurements at larger p t ^ are required. The measured values of v, for which no prediction is 
available from the k t factorisation approach, favour the full NRQCD prediction. 



5 Summary and Conclusions 

An analysis of inelastic photoproduction of J/ip mesons is presented. The kinematic region 
covers 60 < W 1P < 260 GeV, p\^> 1 GeV 2 and 0.05 < z < 0.9. Cross sections in the low 
z region, z £ 0.3, are presented for the first time. The data can be described by boson gluon 
fusion calculations. In the low z range the agreement between data and (LO) calculations is im- 
proved by including resolved photon processes although the uncertainties due to contributions 
from b decays are appreciable in this region. The differential cross sections are compared with 
calculations in three different theoretical frameworks. 

Firstly, in the medium z range, 0.3 < z < 0.9, (double) differential cross sections are compared 
with calculations in the Colour Singlet Model (CSM) in NLO for direct photons. These are 
found to give a good description of the distributions in W 1P , z and pf ^ (tested up to a value of 



16 



60 GeV 2 ). The NLO calculations show a hard p t spectrum which describes the data well 
in contrast to the LO calculation. The distribution of the transverse momentum of the J/ ip is also 
found to be well described in three separate regions of z. The estimated theoretical uncertainties 
due to the uncertainties in the charm mass and a s are much larger than the experimental errors, 
so the data may be used to constrain parameters of the model. 

A second comparison in this medium z range is made with a LO calculation within the CSM, 
allowing the interacting gluons to have transverse momentum (k t factorisation approach). This 
Monte Carlo calculation is also found to give a good description of the data. 

The third comparison in the whole z range (0.05 < z < 0.9) involves LO non-relativistic QCD 
calculations (NRQCD) including colour octet as well as colour singlet contributions. This is 
of importance since the Colour Singlet Model (in LO) does not describe charmonium produc- 
tion in pp collisions and has fundamental theoretical problems in the description of p-wave 
states. The present differential cross sections can be reasonably well described in shape by the 
NRQCD calculations when including direct and resolved contributions. The normalisation de- 
pends on the details of the calculations, in particular on the chosen octet long distance matrix 
elements (LDMEs). The present photoproduction data clearly favour very low values of the 
octet LDMEs. At high z, resummed NRQCD calculations applicable at large p t ^, tend to im- 
prove the agreement with the data. At low z the inclusion of colour octet contributions improves 
the agreement between data and theory. However, decays of b flavoured hadrons are expected 
to play an important role here and are neglected in the present analysis. 

The p\ , distribution is reasonably reproduced by the NRQCD calculation, which however gives 
a slightly softer p^ dependence than that measured in the medium and in the low z ranges. In- 
cluding NLO effects would probably improve the agreement with the measurements as already 
seen in the CSM. 

The polarisation parameters of the J/ip meson have been measured as a function of p t ^ and z. 
Within present experimental and theoretical errors NRQCD, the CSM and the k t factorisation 
approach all fit the data reasonably well. 

With an appropriate choice of parameters, theoretical calculations within the NRQCD approach, 
constrained by results from pp collisions can describe the present photoproduction measure- 
ments. Next- to-leading order corrections are likely to improve the agreement even after colour 
octet contributions are included. The data have the potential to reduce the current large uncer- 
tainties in the phenomenological parameters. This contributes to the development of a unified 
understanding of charmonium production in different environments such as pp, ep and ^p col- 
lisions. 
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Parameters NLO CSM (M. Kramer M) 



Figures |] and |5 
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Renormalisation/Factorisation scale 

(0[l, 3 £i]> 
m c 

a s (M z ) 



MRST 

\/2m c , max[y/2m c , 1/2 + pf\ (forp^) 
1.16 GeV 3 

1.3 < m c < 1.5 GeV 
0.1200 ±0.0025 



Parameters NRQCD (M. Kramer [IZT1]) 
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Resummed calculations (Beneke et al. [|TT]]) 
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Figure [7] 
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(1.3 ± 0.09) GeV 3 
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Table 7: Parameters used in the QCD calculations which are compared to the data. The expres- 
sions "(0[8, 3 Si\) " etc. denote the long distance matrix elements. The colour octet LDMEs are 
extracted from high p t J/ ip production in pp collisions for the NRQCD calculations o//|7^,|27]/ 
and from data on B decays in /[77|7- For more details and for a discussion of the uncertainties 
see these references. 
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Figure 3: Comparison between data and Monte Carlo simulations after all selection cuts and 
after subtraction of non-resonant background, a) W w , b) z and c) p\^ distributions in the 
range 0.3 < z < 0.9, 60 < W 1V < 240 GeV and pjL > 1 GeV 2 (dataset I), d) W 1P , e) z and 
f) pjL distributions in the range 0.05 < z < 0.45, 120 < W 1P < 260 GeV and pf ^ > 1 GeV 2 
(dataset II). The error bars on the data points are statistical only. The EPJPSI simulation of 
direct + resolved photon processes (full histograms) and resolved photon processes (shaded 
histograms) are compared with the data (for the relative normalisation see text). An estimate of 
the direct photon contribution from b — ► J /if) + X from EPJPSI ( which is not included in the 
full histogram) is also shown as a dotted histogram in d)-f). 
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Figure 4: Inelastic J/ip production in the region 60 < W 1P < 240 GeV, 0.3 < z < 0.9 and 
p\^> 1 GeV 2 . a) Total cross section as a function ofW lv , b) differential cross sections da/dz 
and c) da /dpi ^. The HI data are shown together with NLO calculations in the Colour Singlet 
Model /|7?|/ with MRST /p5|7 parton density functions (W 7P = 115 GeV for b) and c)). The 
shaded band reflects the uncertainties in m c and a s (see Table 0); the dashed (dash-dotted) 
curve is calculated with m c = 1.3(1.4) GeV, a s {M z ) = 0.1175(0.1225). In c) a CSM LO 
calculation with MRST parton distributions, m c = 1.3 GeV, a s (Mz) = 0.1225 is also shown 
(dotted). 
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Figure 5: Double differential cross section dV/ dpf ^ dz m ?/ze same kinematic region and 
pared with the same CSM NLO calculations as in Fig. 
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Figure 6: Differential cross section da / dpi ^ for inelastic J/ip production in two kinematic 
regions, a) 0.05 < z < 0.45, 120 < W 1P < 260 GeV; b) 0.3 < z < 0.9, 60 < W JP < 240 GeV 
(same data as in Fig. f^c). The data are compared with LO theoretical calculations in the 
NRQCD framework, including colour octet and colour singlet contributions from direct and 
resolved photons shown as the shaded band, which reflects the uncertainties due to the LDMEs. 
The dashed line shows the colour singlet contribution separately. In b) the band of the same 
CSM NLO calculation /|7?|/ as in Fig. f?]c w also shown as dash-dotted lines. 
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Figure 7: Differential cross section da/dzfor 120 < W rp < 260 GeV and pf ^ > 1 GeV. The 
data are shown as triangles (points) corresponding to datasets II (III). The two data points from 
different analyses at z « 0.4 are statistically correlated and offset for visibility. Theoretical cal- 
culations within the NRQCD/factorisation framework including colour octet and colour singlet 
contributions from direct and resolved photons are shown for comparison. The shaded band is 
a calculation by Kramer /|27J/ reflecting the uncertainties due to the LDMEs ( see Table 0). The 
two dash-dotted lines indicate the region of the resolved contributions (CS+CO) separately and 
the dotted line shows the colour singlet contribution. The dashed line is the result of a NRQCD 
calculation by Kniehl et al. /|7^/, where higher order effects have been estimated and which has 
been normalised to the data. Note that the theoretical calculations are for charm only, while 
the data contain a residual background from b quarks at low z. 
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Figure 8: Differential cross sections da/dz (60 < W 1P < 240 GeV)for a) p t ,ip > 2 GeV and 
b) Pt,ip > 3 GeV in comparison with NRQCD calculations including colour octet and colour 
singlet contributions after resumming soft contributions at high z /|77[/. The curves correspond 
to three values of the parameter A = (no resummation), 300 and 500 MeV (dashed, full and 
dash-dotted lines). The dotted line is the colour singlet contribution alone. The theoretical 
curves have been scaled with a common factor Una) and 3 in b), respectively. 
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Figure 9: Inelastic J /ip production in the region 60 < W 1P < 240 GeV, 0.3 < z < 0.9 and 
pl^> 1 GeV 2 in comparison with a k t factorisation model implemented in the Monte Carlo 
generator CASCADE /|73]/. a) Differential cross section da/dz, b) a 1P as a function of W 1P for 
0.3 < z < 0.8 and c) da/dp: 
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Figure 10: Differential cross sections 1 / a da /d cos 9* in the rest system of the J/ip meson for 
different z regions (kinematic region 60 < W < 240 GeV, pf ^ > 1 GeV 2 ) normalised for 
| cos#*| < 0.9. The inner error bars indicate the statistical uncertainty, while the outer error 
bars include the statistical and systematic uncertainties added in quadrature. The shaded bands 
reflect fits to the form ~ 1 + A cos 2 9*. The contours o/±l standard deviation are shown. 
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Figure 11: Polarisation parameters X (a and c) and v (b and d) in the target frame as functions 
of z (a and b) and p t ^ (c and d). The error bars on the data points correspond to the total 
experimental error. The theoretical calculations shown are from the NRQCD approach /[7P|/ 
(shaded bands) with colour octet and colour singlet contributions, while the dashed curves 
show the result from the colour singlet contribution separately. For A a calculation within the 
k t factorisation approach /|27j/ is also shown (dotted line). 
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W 1P [GeV] 


(W JP ) 


0.3 < z < 0.9 


[nb] 

0.3 < z < 0.8 


60-80 


69.3 


23.0 ±2.5 ±3.2 


17.5 ±2.0 ±2.4 


80 - 100 


89.4 


29.7±2.2±4.1 


21.4 ± 1.7±3.0 


100- 120 


109.5 


34.1 ±2.3 ±4.7 


27.2 ±2.0 ±3.8 


120- 140 


129.6 


36.7±2.6±5.1 


27.7±2.2±3.9 


140- 160 


149.6 


37.3 ±3.0 ±5.2 


27.0 ±2.4 ±3.8 


160-180 


169.6 


42.0 ±3.9 ±5.8 


31.8 ±3.1 ±4.4 


180-210 


194.2 


41.7 ±4.9 ±5.8 


31.1 ±3.7± 4.3 


210 - 240 


224.2 


+5 9 

37.8 ± 7.8 _to 9 


+4.1 
26.4 ±5.8 \ 
—3.5 



z 


(z) 


p t ^ > 1 GeV 


da 7P /dz [nb] 
Pt ^> 2 GeV 


Pt,$ > 3 GeV 


0.30-0.45 


0.38 


27.5 ±2.4 ±3.8 


11.6 ± 1.4 ± 1.6 


4.3 ±0.8 ±0.6 


0.45-0.60 


0.53 


51.7 ±2.8 ±7.2 


20.8 ±1.6 ±2.9 


7.3 ±0.8 ±1.0 


0.60-0.75 


0.68 


64.3 ±3.4 ±8.9 


26.8 ±2.0 ±3.7 


10.0 ± 1.1 ± 1.4 


0.75-0.90 


0.83 


77.4 ± 5.0 ± 10.8 


28.3 ±2.6 ±3.9 


10.5 ± 1.3 ± 1.5 



PtJGeV 2 } 


<P?,*> 


d<7 7P /dp t % [nb/GcV^] 


1.0 - 2.1 


1.54 


9.45 ± 0.61 ± 1.31 


2.1-3.5 


2.78 


5.34 ±0.37 ±0.74 


3.5-5.4 


4.38 


2.77 ±0.22 ±0.38 


5.4-7.6 


6.43 


1.46 ±0.13 ±0.20 


7.6-10.0 


8.74 


0.887 ±0.088 ±0.123 


10.0-13.5 


11.6 


0.452 ± 0.053 ±0.063 


13.5-20.0 


16.4 


0.160 ± 0.022 ±0.022 


20.0- 26.5 


23.0 


0.0815 ±0.0120 ±0.0113 


26.5 - 40.0 


32.3 


0.0236 ± 0.0051 ±0.0033 


40.0-60.0 


48.5 


0.0051 ±0.0016 ±0.0007 



Table 8: Medium z: Cross sections in bins ofW lp , z and pi ^ (0.3 < z < 0.9, 60 < W 1P < 
240 GeV). 





(pU) 


0.3 < z < 0.6 


d 2 <7 7P /d^cb [nb/GeV z ] 
0.6 < z < 0.75 


0.75 < z < 0.9 


1.0 - 2.1 


1.54 


12.2 ± 1.2 ± 1.7 


14.0 ± 1.7 ± 1.9 


22.6 ±2.8 ±3.1 


2.1-3.5 


2.78 


5.93 ±0.65 ±0.82 


11.9 ± 1.33 ± 1.66 


11.4 ± 1.8 ± 1.6 


3.5-5.4 


4.38 


3.11 ±0.41 ±0.43 


6.14 ±0.77 ±0.85 


6.27 ± 1.03 ±0.87 


5.4-7.6 


6.43 


2.08 ±0.26 ±0.29 


2.61 ±0.46 ±0.36 


2.72 ±0.57± 0.38 


7.6-10.0 


8.74 


1.17 ± 0.18 ±0.16 


1.56 ±0.27 ±0.22 


2.57 ±0.54 ±0.36 


10.0-13.5 


11.6 


0.530 ±0.103 ±0.074 


1.14 ±0.208 ±0.158 


1.01 ±0.24 ±0.14 


13.5 - 20.0 


16.4 


0.184 ± 0.039 ±0.026 


0.547 ±0.113 ±0.076 


0.276 ±0.083 ±0.038 


20.0- 26.5 


23.0 


0.145 ± 0.029 ±0.020 


0.143 ±0.040 ±0.020 


0.110 ±0.039 ±0.015 


26.5 - 40.0 


32.3 


0.0323 ±0.0097 ±0.0045 


0.0296 ±0.0169 ±0.0041 


0.0804 ±0.0284 ±0.0112 


40.0-60.0 


48.5 


0.0113 ± 0.0059 ±0.0016 


0.0059 ± 0.0029 ±0.0008 


0.0057 ± 0.0033 ±0.0008 



Table 9: Medium z: Double differential cross sections in z and pi ^ (60 < W 1P < 240 GeVJ. 
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W 1P [GcV] 


(W JP ) 


<7 7p [nb] 


120 - 150 


133.8 


12.1 ±1.5 ±2.6 


150 - 180 


163.9 


12.8 ±2.1 ±2.5 


180 - 260 


213.7 


15.6 ±2.5 ±3.1 



P t 2 JGeV 2 ] 




da 7P /dp t % [nb/GeV^] 


1.0 - 2.1 


1.54 


3.56 ±0.83 ±0.68 


2.1-3.5 


2.78 


1.85 ±0.42 ±0.35 


3.5 - 5.4 


4.38 


1.24 ±0.30 ±0.24 


5.4-7.6 


6.43 


0.875 ±0.177 ±0.167 


7.6-10.0 


8.75 


0.489 ±0.115 ±0.095 


10.0-13.5 


11.6 


0.267 ± 0.068 ±0.053 


13.5 - 20.0 


16.4 


0.117 ±0.030 ±0.024 


20.0-40.0 


27.8 


0.0197 ±0.0044 ±0.0043 



Table 10: Low z: Differential cross sections in bins ofW lp and p^ (0.05 < z < 0.45, 120 < 
W / 7P < 260 GeVJ. 



z 


(z) 


da^p/dz [nb] 


0.05-0.15 


0.10 


35.1 ±8.3 ±8.5 


0.15-0.30 


0.23 


31.1 ±4.1 ±6.2 


0.30-0.45 


0.38 


36.7 ±4.0 ±7.0 


0.30-0.45 


0.38 


31.8 ±3.6 ±4.4 


0.45-0.60 


0.53 


58.7 ±4.3 ±8.2 


0.60 - 0.75 


0.68 


80.2 ±5.3 ±11.2 


0.75-0.90 


0.83 


98.7 ±7.5 ±13.7 



Table 11: Total z range: Differential cross sections in bins of z (120 < W 1P < 260 GeV and 
p t ,i> > 1 GeV J. The points at z = 0.375 are statistically correlated, the first one is obtained 
from dataset II the second one from dataset III. 
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